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Abstract 

The simplest field theory description of the multivariate statistics of forward 
rate variations over time and maturities, involves a quadratic action containing 
a gradient squared rigidity term. However, this choice leads to a spurious kink 
(infinite curvature) of the normalized correlation function for coinciding maturities. 
Motivated by empirical results, we consider an extended action that contains a 
squared Laplacian term, which describes the bending stiffness of the FRC. With the 
extra ingredient of a 'psychological' future time, describing how the perceived time 
between events depends on the time in the future, our theory accounts extremely 
well for the phenomenology of interest rate dynamics. 

Forward interest rate f{t, x) is the interest rate, agreed upon at time t, for an instan- 
taneous loan to be taken at future time x > t, between x and x + dx. At any instant 
of time t, the forward rate curve (FRC) f{t, x) defines a kind of string which moves and 
deforms with time. Modelling the motion of this curve is of paramount importance for 
many financial applications [T]: pricing of interest rate derivatives, such as 'caps' that 
limit the rate of loans that individuals take on their housing, interest swaps, risk man- 
agement (asset liability management), etc. The industry standard is the so-called HJM 
model [3 EI- This model has recently been generalized in different directions [1101^], in 
particular by one of us [H El EI, who has proposed a two-dimensional Euclidean quan- 
tum field theory for modelling the forward interest rate curve. The forward interest rate 
dynamics has a drift velocity a{t,x) and volatility a{t,x); it is convenient to define a 
driftless noise field A{t, x) by 

df 

-^it,x) = a{t,x) + ait,x)A{t,x). (1) 



The noise field describes the external shocks of the economy on the different maturities of 
the forward rates; its statistics is governed by the exponential of an 'action' S[A], which 



gives the weight of a given path of A in the two dimensional space x, t, and is defined on 
the semi-infinite domain x > t. The simplest action that factors in the one dimensional 
nature of the forward rate string was proposed in [H |Hl El , and reads 
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where /i is a 'rigidity' parameter, that gives an elasticity to the FRC To eliminate bound- 
ary terms from the action we choose to impose Neumann boundary condition, i.e. 



dA{t,x) 



dx 



x=t 



0, (3) 



corresponding to a parallel motion of the FRC for short maturities, which is reasonable 
since the spot rate f{t, t) is fixed by the central bank, and very short maturities carry no 
extra risk. 

The quantum field theory [10] needed to describe the statistics of the FRC is defined 
by a functional integral over all variables A{t, x) and yields the partition function Z = 
J DAe^^'^\ The propagator (or noise correlator) is given by 

{A{t, x)A{t', x')) = ^ J DA A{t, x)A{t', x')e^^^^ = 5{t - t')D{x, x'; t). (4) 

Since the above action is Gaussian, the market defined by the above model is complete in 
the sense that all contingent claims can be perfectly replicated by hedging appropriately, 
as in the standard Black-Scholes model. For consistency of the description, one should 
further impose a 'martingale' condition that reads [7j 



a{t, 9) = a{t, 9) [ d9'D{9, 9')a{t, 9'). 
Jo 



Note however that this term is usually numerically very small . With the above choice 
of the action, the propagator, in new co-ordinates 9± is given by 

D{9+,9.) = ^[e-'^«+ +e-'^l^-l] 

(5) 

with 9± = 9±9', 9 = X -t and9' = x' -t with x, x' > t. 

Note that the slope of the propagator perpendicular to 9_ = 0, as in Figure [TJ is 
discontinuous across the diagonal 



dD(9+;9^ 

m = 



d9_ 



/i' f -1 > 

fl 9_<0 



All the variants of the propagator based on a gradient squared rigidity term in the 
action show a similar infinite curvature singularity along the diagonal [TT] . 

However, as discussed in [TT], the surface of the empirical propagator given in Figure 
121 is extremely smooth and shows no such kinks. This observation is in fact related to the 
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Figure 1: For 6_ = 9 — 6', the figure shows that the diagonal axis is given by 6^ = 0. 
The direction of change in 6_ for constant 6*+, namely 66 is orthogonal to the diagonal, 
as shown in the Figure. 

empirical study of (see also [SI ^3)' where the correlation of the (fixed maturity) 
forward rates daily variations 6f{t,9) = f{t + e,t + e + 9) — f{t,t + 6) was studied. More 
precisely, the eigenvectors '^qiO) and eigenvalues Cg, g = 1, 2, of the correlation matrix 
M{e,e'), defined as 

N{e,e') = <6fit,e)6fit,e')>, 

were determined. It was found that the eigenvectors show a structure similar to the 
modes of a vibrating string (\E'i has no nodes, has one node, etc. - see also ^15j). 
and that the eigenvalues (g behave as (a + hq^)^^ for small q (where a, h are constants), 
crossing over to a faster decay ~ q~^ for larger g's. It is clear that the square gradient 
action does indeed lead to plane- wave eigenvectors, with eigenvalues given by {a + hq^)~^ 
with h oc However, the q~^ falloff points to the existence of another term in the 

action that is a fourth power in the derivative along the future time x. This higher power 
of the derivative stiffens the fiuctuations of the forward rates curve, in the sense that two 
nearby maturities experience more correlated external shocks. The aim of the rest of this 
letter is to compare in details empirical data with the predictions of the following action 

that includes the new stiffness term. We will see that this theory indeed allows one to get 
rid of the infinite curvature of the propagator along the diagonal, but that a quantitative 
agreement with empirical data can only be achieved if the action is written in terms of 
an effective 'psychological' time z in the maturity direction, that is a sublinear function 
of the 'true' maturity x — t. In other words, changes in the maturity direction are given 
by d/dz. The introduction of z is similar (but not equal [TT]) to the rigidity and the 
stiffness A^^ not being constant along the x direction, which should be expected. With 
this extra ingredient, we will see that the details of the propagator surface are reproduced 
with surprising accuracy. 



Observed Correlation 




Fiffure 2: Empirical correlation C(6,6') = {Sf{t,0)Sf{t,e ))e ^^ determined from the 
Eurodollar rates in the period 1994-1996. See pT^j for more details on the data. 

Let us first compute the propagator of the 'stiff' action above. Using eq.dHJ, we find 



q;_i_ — a;_ / a_ a 
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with 



1±a/1-4(5) 
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In the limit A ^ oo one finds ^ -^^//^^ and a_ ~ /i^. Hence the propagator has the 
following limit 

lim G{e+; e^; /i, A) ^ D{e+; 9^; fi) 

A— »oo 

and, as expected, reduces to the case of constant rigidity. 

The solution for a± yields three distinct cases, namely, when a± is real, complex or 
degenerate depending on whether fj. < v^A, /i > V2\, ji = ^/2X respectively. One finds 

__A_ ^g-A0+cosh(fe) gi^j^i^ ^ sinh(fe)} + e-^l^-l™'^'^^ sinh{fe + A|e_| sinh(6)}] 
A [e-^e+{i + X9+} + e-^l^-l{l + X\9^\}] (9) 

k 2sin(2<A) [e-^^+"°'^W sin{0 + Ag+ sin(0)} + e~^l^-l"°'^W sin{0 + A|g„| sin(0)}] 

In the above equation, a± = A^e^'' in the first case, and a± = A^e^*"^ in the third case, 
and 6 = = in the degenerate case. Expanding the propagator G{9+; 6 J) about 6*- = 
leads to a cancellation of the term linear in |6'_| and gives a final result that is a function 
of More precisely, the curvature tq orthogonal to the diagonal line 6'_ = is given, 
in the real case, by 

d^Gb{e+-e^^ A3sinh(6)r ,2/,x w,M / X 



rg < follows from the fact that b > 0, confirming that the value of the propagator 
along 5_ = is a maximum. A similar result holds in the complex case. Note that in 
the limit of A ^ cxd, one can no longer carry out the Taylor expansion around 6^ = 0, 
the cancellation of the term linear in |^_| becomes invalid, and the propagator G{6+; 9-) 
develops the expected kink. 



In order to compare with empirical data, we define the normalized correlation function 



by 



^G{e,e)G{e',e') 



Of special interest will be the curvature of C{6, 9') perpendicular to the diagonal, for which 
an explicit expression can be obtained. We do not show this formula here but note that 
the curvature is predicted to increase with 6*+ (see Fig El inset). As one moves along the 
diagonal to longer maturities, the (negative) curvature of C{9, 9') increases, which means 
that the noise affecting nearby maturities is faster to decorrelate as a function of 9 — 9'. 
This is contrary to intuition: since the long term future is much more uncertain, one feels 
that shocks in the far future are more difficult to resolve temporally than shocks in the 
near future. Therefore one expects, and indeed empirically find, that the curvature is a 
decreasing function of the maturity. In fact, we have discoved the new following 'stylized 
fact': the curvature R of the FRC correlation function along the diagonal decays as a 
power law of the maturity, R{9+) ~ 9^^ ^ with v 1.32 (see FiglHl inset). 




Figure 3: Circles: Empirical correlation C{9] 9') along the longest stretch perpendicular 
to the diagonal, i.e. 9' = Q — 9, where G is the maximum available maturity. The plain 
line is the best fit with our stiff propagator model and a power-law psychological time. 
The inset shows a plot, in log-log coordinates, of the curvature —R{9+) as a function of 
9^, for a) the stiff model in real (physical) time, showing an increasing curvatue (dashed 
line) and b) the stiff model with a power law psychological time z{9), correctly sloping 
downwards as to reproduce the power law behaviour of the empirical curvature (squares). 



A way to reconcile the above finding with our stiff quantum field theory is to realize 



that the gradient terms in the action need not be uniform. More intuitively, in the mind of 
market participants, the perceived time between events depends on the time in the future, 
and a decreasing function of the maturity itself. The distance (in time) between - say - a 
10 years maturity and a 30 years maturity is clearly much less than the distance between 
a 1 month maturity and a 10 year maturity. A way to describe this mathematically is to 
replace the true, physical future time 6 = x — thy a. psychological time z = z{6) [TTl lltij. 
which is expected to grow sublinearly with 6, so that the rate of time change z'{6) is a 
decreasing function of 6. In general, one can impose some general features of function 
z{6): it is monotone increasing, such that 6 = 9{z) is well defined, and one can impose 
z{Qi) = ; z{oo) = oo. The independent variables will now be t,z{9) instead of t,x. Our 
final model for the forward rate dynamics then reads: 

^it,e) = ait,zie)) + ait,zie))Ait,zie)) e = x-t (12) 

where f{t,6) depends only on calendar time 9 = x — t. An important feature of the 
defining equation above is that both future times, namely 9 = x — t and psychological 
time z{9) occur in the theory ^. The corresponding stiff action in psychological time is 
written as 

T rco poo / l//:5/l\2 T /^2a\'^\ 

(13) 

The propagator for is G{z, z'; fi, A) as in eq.® and the martingale condition for psy- 
chological time is given by |7j a{t, z) = a{t, z) ^^^-^ dz'G{z, z')a{t, z'). Our final result on 
the normalized correlation now reads 

Cc,.{9+;9_) = 9Az.) + 9-{^-) ^^^^ 

^[g+iz+ + z^)+ g-{0)][g4z+ - z-) + (7-(0)] 
z±i9+;9^) ^ zi9)±zi9') 

with, in the real case that will be of relevance for fitting the empirical data 

g+{z) = e-^"^°'^'^('')sinh{6 + A2sinh(6)} 
g_{z) = e-^l^l™'^'^(^)sinh{6 + A|2|sinh(6)} 
A2 
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It can be shown [7] that the curvature with psychological time reads 

d'CQ,{9v,0) 
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[z\9^)fRQ{2z{9^/2)), (15) 



where Rq is the curvature of the model in physical time. Since one observes a power law 
fall off for the curvature, we can make the ansatz z{9) = 9^ for fitting the data. Using the 

^The theory for psychological future time can be defined entirely in terms of forward rates f(t, z{9)). 
However, for imposing the martingale condition, it is necessary [7] to specify the relation between 
f{t, z{9)) and /(i, 6), and in effect one would recover Ea. H12l) 



fact that Rq{2z{9^/2)) varies very slowly as a function of 6'+, one can make the following 
approximation 

[z\e^)f^^^^2^-2 = -u, (16) 

leading to 77 ~ 0.34. Therefore the psychological time flows, as expected, much slower 
than real time. The rate of change of psychological time decreases as ~ ^-o-66. jQaj: 
after ten years looks similar to two weeks after a month. 

Having used the behaviour of the curvature to fix the value of rj (and thus, up to 
an irrelevant overall scale, the function z{9), we are left with only two parameters, A 
and /i, to fit the whole correlation surface C{d,6'). For the Eurodollar data that we 
have used (see for more details), we have up to 30 different maturities and therefore 
405 different points (the diagonal values are trivial). We determine A and /i such as to 
minimize the average square error between the empirical C{6, 9') and the prediction of 
the model. Defining A'' = A and jl^ = /i, the best fit is obtained for A = 1.79/year and 
jl = 0.403/year, corresponding to 6 = 0.845. The residual error is as low as 0.4% per 
point, and the order of magnitude of the time scales (one year) defined by A and fl are 
very reasonable. The remarkable quality of the fit can be checked in FigOl along the 
longest stretch perpendicular to the diagonal, i.e. 6' = Q — 6, where O is the maximum 
available maturity (7.5 years). Even more remarkable is that the curvature of C{6,6') 
along the diagonal is very precisely reproduced by the same fit, as shown in the inset of 
Fig. El Testing the fit on the second derivative of the fitted surface is of course much 
more demanding. The existence of the boundary at x = t, or ^ = 0, is refiected in the 6+ 
term in the propagator; if one removes this term, and in effect assumes that the forward 
rates exist for all —00 < x < +00, then the fit deteriorates with the root mean square 
error increasing from 0.40% to 0.53%. The existence of the boundary at x = t can hence 
be seen to have a significant effect on the correlation of the forward rates. 

Let us summarize what we have achieved in this study. The simplest field theoret- 
ical description of the multivariate statistics of forward rate variations over time t and 
maturities 6, involves a quadratic action containing a gradient squared rigidity term jH], 
that captures the one dimensional string nature of the forward rate curve. However, this 
choice leads to a spurious kink (infinite curvature) of the normalized correlation function 
along the diagonal 6 = 9'. Motivated by a previous empirical study [12] where the short 
wavelength fiuctuations of the FRC were shown to be strongly reduced as compared to 
that of an elastic string, we have considered an extended action that contains a squared 
Laplacian term, which describes the bending stiffness of the FRC. In this formulation, 
the infinite curvature singularity is rounded off. In order to fit to the observed correlation 
functions of the forward interest rates, however, one has to add as an extra ingredient 
that the rigidity/stiffness constants are in fact not constant along the maturity axis but 
increase with maturity. An intuitive and parsimonious way to describe this effect is to 
postulate that markets participants, who generate the random evolution of the FRC, do 
not perceive future time in a uniform manner. Rather, time intervals in the long term 
future are shrunk. The introduction of a 'psychological time' z{9), found to be a power 
law of the true time, allows one to provide an excellent fit of empirical data, and in 
particular to reproduce accurately a new stylized fact: the curvature of the forward rate 
correlation perpendicular to the diagonal decays as a power-law of the maturity. We 
believe that our quantum field formulation, including the new stiffness term and coupled 
with an appropriate calibration of the term structure of the volatility (j{9) (see [T^ I13j). 



accounts extremely well for the phenomenology of interest rate dynamics. It is also math- 
ematically tractable, and should allow one to compute in closed forms derivative prices 
and optimal hedging strategies. It would be interesting to generalize the above model 
to account for non Gaussian effects, that are important in many cases fl]. This would 
amount to considering non quadratic terms in A in the action. We leave this extension 
for future work. 
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